Transverse quark motion inside charmonia in diffractive photo- and electroproductions 
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We reexamine the Fermi motion effects on the diffractive photo- and electroproductions of heavy 
vector-mesons, J/tp and ip' j oS a nucleon in the leading ln(Q^/AQCD) approximation (LLA) of 
QCD. We take into account all the Fermi motion corrections arising from the relative motion of 
quarks inside the charmonium, which is treated as a nonrelativistic bound state of c and c. Our 
key ingredients are the correct spin structure for the cc bound state, the off-shellness in the 
cc — » J/tp {tp') hadronization vertex, and the modification of the gluon's longitudinal momentum 
fraction probed by the process, due to the relative motion between c and c. We demonstrate that 
these three contributions produce the new Fermi motion effects in the LLA diffractive amplitude in 
QCD. It is found that our new effects moderate the strong suppression of the diffractive J/tp (tp') 
production cross sections, which was reported in the previous works on the Fermi motion effects. We 
emphasize the role of the transverse quark motion for the heavy meson production, and also discuss 
the strong helicity dependence of the Fermi motion effects and its implication in the longitudinal to 
transverse production ratio, ol/ot- 
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I. INTRODUCTION 

High-energy diffractive photo- and electroproductions 
of charmonia off a nucleon, 7*^*^ + N^V + N(y^ 
J/tp, tp'), are of particular interest, since they open a new 
window on the interface between perturbative and non- 
perturbative aspects in QCD With perturba- 

tive QCD (pQCD), it has been shown |2| that the cross 
sections of these processes depend quadratically on the 
relevant hadronic matrix elements, i.e., the wave function 
(WF) of a charmonium {J/tp or tp') and the gluon dis- 
tribution inside the target nucleon. Thus, the diffractive 
charmonium productions are expected to be sensitive to 
those nonperturbative matrix elements, especially the be- 
havior of the gluon distribution for small Bjorken-x. Re- 
cent data from the ep collider HERA 5, 6, 7, 8, Q.ITollTll 
will indeed provide us with quantitative information on 
them, while the precision of the data calls for QCD anal- 
ysis taking into account also the subleading effects. 

For high 7^*-'-A^ center-of-mass energy W , the corre- 
sponding diffractive amplitude obeys factorization and, 
in the leading logarithmic approximation (LLA) of 



pQCD, is expressed schematically as (see Fig.^ 



A1 = 



ccN 



(1) 



with the photon light-cone WF ^'^ for the 7'*^ — > cc 
transition, the amplitude A^'^^ for the elastic scattering 
of the cc-pair off the nucleon by exchanging reggeized 
two gluons, and the vector-meson WF for the soft 
hadronization process cc V 12. Here, A"""^ is further 
factorized into the cc-gluon hard scattering amplitude 



jccg -(-j-^g gluon distribution G inside the nucleon. 

Physically, the factorization of Eqs. and follows 
from the property that the typical timescale of the cc-iV 
scattering at high W is much shorter than the lifetime of 
the cc fluctuation, as well as of the fluctuation of partons 
inside the nucleon . Moreover, the heavy-quark mass 
nic ensures that the 7'*-' — > cc fluctuation is sufficiently 
a short-distance process that takes place within a dis- 
tance of order X/rxic or less. This enables us to treat the 
photon WF , as well as the scattering amplitude A'^'^^ 
between the small-size cc-color-dipole and the gluon, by 
perturbation theory even for the photoproduction with 
the photon's virtuality = 0. 
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FIG. 1: Schematic diagram for the diffractive photo- or elec- 
troproduction of a charmonium V{= J/tp,tp') off a nucleon 
iV, 7'*) +N + N. 

Eq. involves the convolution with the cc-bound- 
statc WF Here, the large mass of TOc also allows us 
to organize the amplitude M into a hierarchy in terms 
of the dependence on the average velocity v of the heavy 
quark in the charmonium rest frame. The leading term 
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independent of v corresponds to the static or nonrelativis- 
tic limit, where c and c have no relative motion inside 
the charmonium, while the subleading terms in powers 
of V express the Fermi motion corrections and relativistic 
corrections. For the charmonium, quark potential model 
calculations as well as lattice QCD simulations 
indicate that v is reasonably small, ~ 0.2 ^ 0.3. Thus, 
the subleading terms are not expected to produce so large 
corrections, though they could be physically important. 

So far, there exist several investigations for the Fermi 
motion corrections, but significant discrepancies are ob- 
served among the previous results: Frankfurt-Koepf- 
Strikman (FKS) Q predicted quite strong suppression 
of the cross section for the diffractive photo- (electro- 
)production of J/-0 by a factor < 0.1 (< 0.3) due to the 
Fermi motion effects. On the other hand, the results of 
Ryskin-Roberts-Martin-Levin (RRML) 0| gave a signifi- 
cantly weaker suppression factor of 0.4 ~ 0.6 for the Jjjb 
photoproduction. Furthermore, the author of Ref. Il4|| 
estimated an even weaker suppression factor of 0.9 or 
more for the J/tp electroproduction. The treatment of 
these three works is different in several points; in par- 
ticular, the difference in the detailed shape of the meson 
WF 4"^ appears to be an important source of the dis- 
crepancies, as claimed by the authors of Ref. 01 . In our 
opinion, however, any of those works is not satisfactory 
as the calculations of the Fermi motion effects, because 
certain effects due to the relative motion of the quarks 
are not taken into account. 

The aim of this paper is to reexamine the Fermi motion 
corrections on the diffractive photo- and electroproduc- 
tions of charmonia. Similarly to the previous works 
our basic formulation relies on the QCD factorization ap- 
proach in the LLA of pQCD. In this framework, we re- 
derive the formulae Q explicitly, paying attention 
to the contributions that are related to the relative mo- 
tion of the quarks and thus should be regarded as a part 
of the Fermi motion effects. We find three types of such 
additional contributions, which have not been considered 
properly in a single treatment so far. 

First of all, from the viewpoint of the systematic ex- 
pansion of the total amplitude (^3) in terms of the velocity 
V, the Fermi motion corrections through the relative or- 
der v"^ should be computed with the Bethe-Salpeter cc 
WF VP^, which is defined as the matrix element of a 
bilocal operator composed of the two-component Pauli 
spinor fields for c and c in a nonrelativistic Schrodinger 
field theory, replacing the original relativistic field the- 
ory. Namely, when the WF is expressed in the usual 
four-component Dirac representation, it is subject to a 
constraint TZ'i'^ = with the projector TZ onto the 
"large components" corresponding to the Pauli spinor. 
Physically, this implies that the WF describes a pure 
•^S*! cc-state. RRML did not consider the role of such 
projector TZ at all In FKS 3], an improper projector 
TZ was used, which was constructed in an oversimplified 
fashion. We employ the projector TZ, which obeys spin- 
symmetry in the Schrodinger field theory and coincides 



with that appeared in Ref. 0]. We will show that the 
spin structure of ^E"^, which is imposed by our TZ, modi- 
fies the amplitude ^ at the order in comparison with 
RRML or FKS. 

Following the previous works when combining 

the factors of Eq. we express the meson WF ^E"^, as 
well as the other parts A'^'^'^ and '^'^ , in a helicity repre- 
sentation. This corresponds to a familiar representation 
of the QCD factorization formulae in the Lepage-Brodsky 
approach for hard exclusive processes jl5j. The helicity 
representation of v]/^* in Eq. (^) is given by the matrix 
element v\i'^^*u\, in terms of helicity Dirac spinors ux 
and vx' for quark and antiquark. This spinor matrix el- 
ement representing the effective "cc —^ V vertex" may 
be off the energy shell: the total energy carried by the 
two spinors v\> and u\ is generally different from the en- 
ergy of the vector meson. In the previous works 0, Q , 
however, the on-shellness was assumed when evaluating 
the matrix element for the effective cc V vertex. Ap- 
parently, this assumption corresponds to neglecting the 
binding energy of the charmonium, so that it could give 
the 0{v^) error in the final result. Thus, the second of 
the additional Fermi motion contribution is generated by 
the off-shellness of the cc ^ V vertex. 

Third, the Feynman diagrams relevant for the LLA of 
the factorization formula 10) tell us that the gluon's mo- 
mentum fraction inside the nucleon, which is probed by 
the diffractive production, depends on the longitudinal 
as well as transverse relative motion of c and c through 
the momentum exchange at the quark-gluon coupling. 
Namely, the Bjorken-a; of the distribution G of Eq. Q 
depends on the motion of the quarks, which produces 
the 0{v^) effects to the amplitude (Q, though the pre- 
vious works 0, 0, seem to have simply neglected this 
contribution. 

As we will demonstrate, the modification of the ampli- 
tude due to the above three types of contributions is 
sensitive to the transverse relative motion of the quarks 
inside charmonia; i.e., the transverse quark motion is no 
less important than the longitudinal one for the produc- 
tion of nonrelativistic bound-state with < rn^. This 
is in contrast with the light-cone dominated processes 
such as the hard diffractive productions of the p me- 
son {Q^ ^ ^qcd)' where the longitudinal quark motion 
plays a dominant role compared with the transverse one, 
and the effects of the transverse quark motioiigive only 
small corrections of higher twist ^ ^qcb/Q^ H- 

By including those three types of additional contribu- 
tions, we will give quantitative estimates of the Fermi 
motion effects. We employ the "shape" of the meson 
WF 'i'^ , which is obtained by solving the Schrodinger 
equation with a realistic potential between the quark 
and the antiquark. We compute the cross sections for 
the diffractive photo- and electroproductions of J/tpi and 
find that the Fermi motion effects give moderate suppres- 
sion. The net effect of our Fermi motion contributions 
is, roughly speaking, similar to that of RRML (j, but 
our results show some new features in comparison with 
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those of RRML as well as of Refs. 0, To demon- 
strate them, we calculate the ratio of the cross sections 
with longitudinally and transversely polarized photons, 
ctl/ctt- Because the contribution of the gluon distribu- 
tion G proves to mostly cancel in this ratio, this quantity 
is sensitive to the first two of the above-mentioned addi- 
tional Fermi motion effects, which are strongly helicity- 
dependent. Another interesting quantity is the ratio of 
Vj' to J/tp cross sections, by extending our calculation to 
the case of a radial excitation tp' . We observe that the 
Fermi motion effects produce characteristic behavior for 
the ratio, reflecting the different shape of the WFs of J /'ip 
and tp' . 

The paper is organized as follows: in Sec. 11, we derive 
the formulae and (jSJ in the LLA with our additional 
Fermi motion effects. After recalling the basic steps to 
derive the factorization formula for the diffractive pro- 
duction, the three types of additional Fermi motion ef- 
fects are explained in detail in Sees. II A-C. The final 
form of our production cross sections is given in Sec. II D. 
Sec. Ill contains our numerical results. In Sec. Ill A, 
we present gl, cft and olIot for the J/V' productions, 
demonstrating roles of each of our three corrections. We 
show the total cross sections for the J jtp photo- and elec- 
troproductions in Sec. Ill B, making a comparison with 
the HERA data, and clarify the connection with the pre- 
vious works. Similar program is carried out for the ?/;' 
productions in Sec. Ill C, and also the implication of our 
results in recent debate |2a . l23ll2^ on the V'' to J/V' ratio 
is discussed. In Sec. IV, a summary and discussion are 
presented. 



II. FORMULATION 

A. The forward diffractive amplitude in QCD 

Let us first describe a basic formulation to evaluate 
the diffractive photo- or electroproduction of the heavy 
quarkonium. We consider the forward diffractive pro- 
duction of the heavy vector-meson, 7'-*^(q) -I- -/V(p) 
V((l + A) + — A), where q and p denote the four mo- 
menta of the photon and the incoming nucleon, respec- 
tively, and A the momentum transfer in the t-channel. 
The tot al cente r-of-mass energy of the j^*^-N system, 
W = \/ [p + qY, is assumed to be much larger than the 
photon's virtuality and the heavy-quark mass. 



> ml, 



(3) 



with 



and ml > A^^D 



We also suppose ^ ~t with t = A^, 
Under such kinematical condition, 
the relevant diagrams for the scattering amplitude in the 
LLA are shown in Fig. [21 The heavy-quark mass mc al- 
lows us to calculate the creation of the cc pair, ^''*'^ cc, 
as well as its time development by perturbation theory 
before the nonperturbative effects to form the quarko- 
nium state V set in. As is well-known, the corresponding 
amplitude is dominated by its imaginary part, which we 




N(p-A) 



(b) 



k-q + l k-q k- q + I- A 




(C) 



k—q + l 



k— q 



k+l 











FIG. 2: Feynman diagrams for the diffractive photo- or elec- 
troproduction of V = J/tp,tp'. The dashed line means the 
discontinuity of the amphtude. 



shall calculate explicitly; the contributions from the real 
part can be incorporated perturbatively. For simplicity, 
we consider the forward limit of the diffractive process, 
t ^ 0, in detail. 

To proceed, we introduce two null vectors g' and p' 
and q' ■ p' = s/2 as 



with q'"^ = p'^ 



q = q p , 

s 



p 



Ml 



q'+p' 



where is the nucleon mass, p-^ 
Q^. In the forward limit t - 



Ml 



- Ml , and s 
0, we also have 



VF2 



Ml 



P 

with Mv the mass of the heavy vector-meson, so that [p— 
A)2 = Ml and [q + A)^ = M^. In Fig. El the momenta 
k for the internal quark lines and I for the exchanged 
gluons have the Sudakov decomposition with respect to 
q' and p': 

aq' + f3p' + k^^, (4) 
W + ^P' + 1±- (5) 



k 
I 



Here, a expresses the fraction of the photon momen- 
tum carried by the charm quark along the "g' direc- 
tion" , while ^ the fraction of the nucleon momentum 
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carried by the gluon along the direction" . We choose 
a Lorentz frame as g"*" = q'^ and p~ — p'~ , so that 
k = {aq^ , j3p^ ,k_\_), k\ = —k\, and similarly for I. 

In the calculation of the discontinuity of the relevant 
diagrams, the on-shell conditions for the quark lines, 
which cross the "cut" as shown in Fig.|2Ia), determine (3 
and ^ as 



/3 



1 



.(1- 



[(l-a)02 + fc2 



1 



(6) 



(7) 



while the integrations over a, fc_L, k, and Zj_ has to be 
performed. Now, formally, the factorization formula fol- 
lows from the observation that the dominant integration 
region in the K-variable is the one of k = 0{l\/s) for 
fixed l±_ in the high-energy limit (PJ, and, in this regime 
f' ~ —l\ <C s, the exchanged soft gluons couple to 
the quarks with an eikonal vertex g^' corresponding to 
a single gluon polarization. The use of the eikonal ver- 
tex immediately leads to the factorization of the photon 
WF '^'^ and the vector-meson WF from liaM. (see 
Eq. 1^). Also, the K-integration reduces to that for the 
nucleon matrix element by Taylor expanding explicit k- 
dependence of the other parts in the integrand, and this 
accomplishes the factorization of the gluon distribution 
along with appropriate decoupling of the color indices; 
this is so-called "A;_L-factorization" |^^, because the in- 
dependence in the integrand has to be retained in the 
present regime. As a result, the unintegrated gluon den- 
sity distribution f{£,,l'±), depending on the transverse 
momentum l± as well as the longitudinal momentum 
fraction ^, describes the nonperturbative part due to the 
nucleon matrix element (Actually, the unintegrated 
gluon density f{S,,l']^) for P ~ involves all multiple 
gluon exchanges, which give the "ladder structure" de- 
scribed by the BFKL evolution, but in the LLA we can 
eventually eliminate / in terms of the conventional gluon 
distribution G as in Eq. For example, we obtain as 
the discontinuity of Fig. OJa) , 

TT^as da f d^kj_d^lj_ f{^,ll) 



/m;. sJo a'' J (27r)4 ^4 
X > ^l^y^^*{a,k^)[uy„{k)/uy,(k + l)] 



E 

AA'A"A'' 



X [^lA"(fc + 0/uA'(fc)]*lfP(a,fc±), 



(8) 



up to the terms suppressed for high s. Here, ag — 
g"^ /{Ait), is the number of colors, and Q and C,' de- 
note the helicities of the photon and the vector meson, 
respectively. uy{k) is the on-shell spinor for the quark 
with the helicity A'. 

The light-cone WF ^/jj^ {a,k±) describes the proba- 
bility amplitude to find the photon of helicity C in a state 
with the minimal number of constituents c and c, which 
carry the longitudinal momentum fractions a and I — a, 



the transverse momenta k± and —k±, and the helicity 
A' and A, respectively. It is given as the lowest order 
amplitude for the photon dissociation to the cc-pair (see 
Fig. 121) , and the corresponding formula is well-known 0] : 



n^i\a,k^) 



^a{l-a)uy{k)i'^ihx{q~k) 
a[l - a)Q'^ + k\ + ml 



(9) 



where v\{q — k) is the on-shell helicity spinor for the 
antiquark. e is the charge of the proton and Cc — 2/3 
gives the charge of the c-quark. The polarization vectors 
e^'') of the photon are defined as 



,(0) 



,(±1) 



Q 

.(±1) 



— p 



(10) 

(11) 



for the longitudinal (C = 0) and transverse {Q = ±1) 
polarizations, respectively. Calculating the spinor matrix 
elements in Eq. with a helicity basis by Lepage and 
Brodsky [l^ , one finds for the longitudinal photon. 



2a(l - a)Q 



.a(l-a)Q2 + fcj 



and, for the transverse photon. 



1 

Q. 



,(12) 



k\+ml + a{l - a)Q'^ 

2eec (e^^ ■ k±){a - Ja.t) 
k^j^+ml + a{l - 0)0^ 



(A' = A) 



(13) 



(A' = -A) 



with T = ±1 and ej^' — "^(Ij^*)- Note that the sec- 
ond term in Eq. H12fl is independent of both a and fcj^; 
as is well-known, this term can be omitted due to cross- 
ing symmetry, since the contribution from this term is 
cancelled by the corresponding contribution from the di- 
agram (c) of Fig. 121 (see Eqs. lO below). 

In Eq. ©, the formation process of the vector meson 
V from the outgoing cc-pair is factorized as the vector- 

meson light-cone WF ^;^///\ {a,k±). Its physical inter- 
pretation is similar to the photon WF \E'j[fP(Q!, fc^), but 



viC) 

A"'A 



{a,kj_) is more sophisticated because it contains 
nonperturbative dynamics between c and c. Its explicit 
formula will be discussed in Sec. II B. 

We note that the eikonal vertex ^' does not change the 
helicity of quark or antiquark for the Lepage-Brodsky 
spinors: 



UAi(fci)/wA2(fc2) 



vxi{ki)i>'vx^(k2 
s 



^1^2 '5ai,A2 



(14) 
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Substituting this into Eq. ((SJ, we immediately get the s- 
channcl hchcity conservation, i.e., Eq. ||SJ) is nonzero only 
for C = C- 

The contributions from the diagrams (b) and (c) of 
Fig. 13 can be computed in a similar way to the dia- 
gram (a). The results are expressed with the same light- 
cone WFs ^^"^ gluon distribution / as in 
Eq. and satisfy the s-channel helicity conservation. 
Combining those contributions with Eq. ©, we get the 
discontinuity of the total diffractive amplitude for the 
helicity C of the incoming photon as 



ImM 



(C) _ 



3s 



dacPk\ 
167r3 



where a'^ 



cN 



x^lf\a,k^)a^'^^l^i\a,k^), (15) 
expresses the operator acting on 



~ccN 



z1 



1 + 1 



2e'^^-^ 1 , (16) 



with b = id/dk±. Three terms in the integrand corre- 
spond to the diagrams (a), (b), and (c) of Fig. [21 respec- 
tively. Note that Eq. (|15|l obeys the standard representa- 
tion of the factorization formula in the Lepage-Brodsky 
approach for the exclusive processes |15ll] . 

In Eqs. ((T3|l . the l± integral converges for 
oo. Moreover, in the integrand, we have the contributions 
behaving like 



1 



1 



Q 



■2 ' 



(17) 



^2 



with Q = mf. + a{l — a)Q , from the denominator of 
^^i-p(a, A;_l) of Eq. Then, the important domain in 



^2 



the ij_ -integral is ij^ ^ Q , and, in particular, a logarith- 
mic contribution ^ ln((5 /■'^qcd) comes from l\ <^ Q . 
Here, up to non-logarithmic corrections, Eq. H16|l can be 

simplified by the replacement 1 — e*'-^ l\b /4, re- 
taining only the leading nonzero term in the power series 
of Z_L • b, which corresponds to the "color-dipole picture" 
[^|2^. In the same accuracy, we neglect the l± depen- 
dence of ^ in f{^,l±) (see Eq. Q), i.e., replace ^ by 



1 

X = — 

s 



fe 



a{l — a) 



Now, using the relation 



2 

xG{x, Q ) 



dV 



- I' 



(18) 



(19) 



with the conventional gluon distribution G, we get, in 
the leading logarithmic accuracy. 



-ccN 



^a4Qls)xG{x,Qle)b\ (20) 



with Q'^g = (Q^ -I- My)/4. This operator corresponds 
to the color-dipole cross section, representing the high- 
energy interaction of a small-size quark-antiquark con- 
figuration with a target. The results (|15|l and l|2U|) 
demonstrate the factorization and (0); note that, to 
derive those results, we have not employed any other 
assumption than the high-energy limit and the LLA. 
Eqs. (|15|) . (Pn|l coincide with the corresponding results 
of Refs. 0,13/3 J except that x of Eq. lfTS|l appears in- 
stead of the usual x = [Q^ -\- My') /s, and that we shall 
use a more precise form for the heavy vector-meson WF 
"^^I'P {a,k±). In the following subsections, we will ex- 
plain that those differences produce the "new" Fermi mo- 
tion corrections of order that were mentioned in Sec. I, 
when expanded in powers of the heavy-quark velocity v 
inside the vector meson. 



B. Spin structure of heavy-meson WF 

The light-cone WF of the vector mesons is defined as 
the Fourier transform of the Bcthe-Salpeter cc WF at 
equal light-cone time x+ — 0. It is a nonperturbative ob- 
ject, which is not fully understood at present. In the same 
spirit as in Refs. [3 01 j we shall construct the shape of 
the heavy- meson WF by solving the Schrodinger equation 
with a realistic potential between quark and antiquark. 
As we will demonstrate below, this indeed gives a conve- 
nient framework to evaluate the Fermi motion effects of 
order v'^ due to the relative motion of the quarks inside 
the nonrelativistic bound state. 

A field-theoretic definition of the corresponding WF is 
provided by the heavy-quark limit of the original Bethe- 
Salpeter WF, i.e., by the Bethe-Salpeter cc WF in terms 
of the two-component Pauli-spinor fields in a nonrela- 
tivistic Schrodinger field theory. Spin symmetry tells us 
that the S'-wave cc states form the quartet consisting 
of the spin singlet and the spin triplet, and the spin 
triplet states describe the vector mesons like J/ip,ip'- 
As is well-known |l7j| . the corresponding coupling of the 
heavy-quark spins is represented by the matrices, when 
expressed in the original four-component Dirac represen- 
tation, 



(21) 



for the spin singlet and triplet states, respectively, with 
the projector TZ onto the "large components" correspond- 
ing to the Pauli spinor: 



n : 



1+f 



(22) 



Here, V denotes the four velocity of the relevant mesons, 
and e*^'»^ denotes the polarization vector of the vector 
meson. Thanks to spin symmetry, the matrices (|21|) 
completely determine the spin structure of the light- 
cone WFs for the pseudoscalar and vector mesons in the 
heavy-quark limit. In a recent work psj . the diffractive rjc 
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and r^f, productions induced by neutrino beam have been 
discussed using the first of the spin WFs |(2U. Now, 
the second one with V = {q + A) /My gives our hght- 
cone WF for the vector meson. In the hclicity (A',A)- 
representation, it reads 

V 1 — Q! My 

(23) 

with the Lepage-Brodsky spinors u\i{k) and v\{q — k). 
The polarization vectors ej^'' are given as 

1 



,(0) _ 



=(±1) 



My 
(±1) 

el 



q p 

s 



(24) 
(25) 



for the longitudinal (C — 0) and transverse (^ = ±1) 
polarizations, respectively, so that e'''' • V = and 
g(C) . ) — — 6(^,(^1. The scalar function (j)*{a,k±) rep- 
resents the nonperturbative part provided by a S'-wave 
solution of the Schrodinger equation, whose detailed form 
will be specified later. We note that the spin structure 
of Eq. reduces to that of the light-cone WF for the 
light vector-meson of Ref. by the replacement TZ ^ I. 

It is straightforward to compute the spinor matrix el- 
ements of Eq. H23|) . which give the effective "cc 
vertex" : 



ri9;(a,fc^) 



1 vxiq-k) 



-r 



ux'{k) 



Mv 

We get, for the longitudinal polarization, 
A(2a-1)(£^^^ -fc^L) 



r 



(0)* 
AA 



r 



(0)* _ 
-AA ^ 



V2 a(l - a)Mv ' 

1 k\ + mcirric + My) 

2 2a(l - a)M^ ' 



V 

(26) 

(27) 
(28) 



where e'^'' = '^(^^ 8i,nd, for the transverse polariza- 
tions (t = ±1), 



Ia* = [Sr,x{{Mv+m,)m, + a{l-a)Mlr} 



V2 a{l - a)Ml 



(29) 



-I- (1 - 2a)mc 



je^l^* ■ fcx) 



(30) 



Using these results, the discontinuity of the total diffrac- 
tive amplitude H15|l can be expressed in the LLA as 



lmM^'^'> ^ J da J d^k±(l)*{a,k±)n^<'>{a,kj 



), (31) 



where fl^^^ (a, k±) represents the "effective 'jNV vertex" , 
whose (a, fc^)-dependence is completely determined by 
the analytic formula: 



J7(«)(a,fe_L) 



3s 



327r3 



^rl9;(a,fc^)a-^vI;5;(C)(o,fc^) 

A'A 



(32) 



with Eqs. dl), ini) for ^1)1' {a,kj_), and Eq. (gOl for 
^ccN ^ Eq. (p?T|l provides the basis of our study for 
the Fermi motion effects: it is given by the integral 
of n'^^^a,k±), weighted by the complicated function 
4>*{a,k±), whose (a, fc_L)-dependence represents the rel- 
ative motion between c and c due to the nonperturbative 
dynamics inside the vector meson. 



C. Fermi motion effects in the LLA diffractive 
amplitude and comparison with other treatments 

For a nonrelativistic cc bound-state, the WF (j>*{a, k±) 
is sharply peaked at a = 1/2 and fc_L = 0. Thus, we may 
conveniently identify the average velocity v of the charm 
quark in the charmonium rest frame as, e.g.. 



V^)c.l2 




with 



_ J da(Pk± ■ ■ ■ (f)*{a, k±) 
J da(Pk±4'* {a, k±) 



(33) 



(34) 



In order to see the dependence of ImA^('») on the veloc- 
ity V, we may Taylor expand il,^'''^{a, k±) around a — 1/2 
and fc^ = in the integrand of Eq. (|31|l . and then per- 
form the integrations over a and k± with the weight 
(f>*{a,k±): the leading term gives the static limit, which 
corresponds to no relative motion of the quarks, i.e., the 
replacement (j)*{a,kj_) — > const x 6{a — l/2)S'^^^{k±^) in 
Eq. (|31|l . and the next-to-leading term gives the Fermi 
motion effects of order w^. Since all realistic nonpertur- 
bative models for 0*(a, kj_) would give similar values for 
Eq. (|33|l as ~ 0.2 ^ 0.3, the size of the Fermi motion 
effects in ImJA^''^ should be essentially determined by 
that of the expansion coefficients in the Taylor expansion 
of fl^''\a, kj_). Namely, the dependence of n''''\a, k±) on 
a and k± determines the actual significance of the Fermi 
motion corrections. 

The (a, fc_L)-dependence of Eq. from ^f^^^p (a, k±) 
is unambiguously determined to be Eqs. (|12|) . (|13|l in the 
LLA. On the other hand, for the other parts r^'5^*(a, k±) 
and a'^'^'^ , the previous works by RRML and FKS em- 
ployed different shapes from our Eqs. H27|) - (|30|l and 
Eq. (pn|l . Because those different shapes could be a ma- 
jor source to influence the Fermi motion corrections, it 
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is instructive to compare our r^'5^*(a, k±) and 0-'='=^ with 
those of the previous works. 

Now, we discuss about the effective cc ^ V vertex 
T^yxi"^ ^-l). In the RRML paper [3|, their effective ver- 
tex reads 



p(0)* _OA 

^ A'A;RRML ~ ^0A',-A 

for the longitudinal polarization, and 
\/2 uIcXStx 



(35) 



^ A'A:RRML ^ ^ 



a{l — a)Mv 



(A' = A) 



(36) 



(A' = -A) 



a{l - a)Mv 

for the transverse polarization with r = ±1. To show 
the relation of the RRML vertex with ours, we calculate 
Eq. H26|l with the replacement TZ ^ 1: 



p(0)* 
^ A'A 



r(±i)^ 

^ A'A 



TZ- 



1 



fc^ -I- ml 



_ p(±i)* 

~ ^A'A;RRML- 



'^A'^-A, (37) 

(38) 



Thus, r;,,;^.i^i^f^L is given by 
the condition, 



.(C)* 

A'A 



with imposing 



+ 



kl 



1 



a 



(39) 



which means that the corresponding cc ^ V vertex is 
on the energy shell: we note that the Lepage-Brodsky 
spinors u\{k), v\{k) are the functions of only the com- 
ponents {k~^,kj_) of the momentum k so that the 
vertex defined as Eq. H26|l obeys the conservation for the 
"three-momentum" {k^,k±) only. Namely, the sum of 
the three-momenta of quark and antiquark equals the 
three-momentum of the vector meson, (g"*" -|- , -|- 
= (9^1 Ox), but the conservation of the "minus" 
component of the momentum, expressed as Eq. (|39|l . 
needs not be satisfied. (See also Ref. 18]. In fact, the sit- 
uation is similar to that for the spinor matrix elements 
of the photon WF lO, that represent the 7 ^ cc ver- 
tex.) Comparing Eqs. (|23J), with Eqs. (123, ijSHI), we 
see that the Fermi motion effects of order would be 
modified if the condition (|39|) were imposed by hand. 

Also, comparing our vertex (|77|l - llSn|) with Eqs. (|S7|l . 
(|38|l . we find that the replacement TZ ^ 1 would modify 
the Fermi motion effects of the order . Apparently, the 
spin structure of Eqs. (|37|1 . (|38l) would be contaminated 
by that for the Z?-wave cc states, so that these vertices 
cannot be combined with nonrclativistic models for the 
WF 0*(a,fc_L). 

On the other hand, FKS also considered their effective 
cc ^ V vertex on the energy shell with Eq. (|39|l . but 
tried to include the spin structure for a pure ^Si cc-state 
They employed the projector, defined in the helicity 
representation as (see Eq. (I26f) ') 

7^FKS = '5c,o'5a',-a + {Sc,i + Sc.-i) 6y^x, (40) 



so that the FKS vertex reads 



(0)* 

A'A;FKS 

(±1)* 

A'A;FKS 



= r 



(0)* 

A'A;RRML' 
(±1)* 



A'A;RRML- 



(41) 
(42) 



The FKS projector l|in|l would give the exact spin struc- 
ture for the triplet states, if it acted on the helicities of 
the two-component Pauli spinors. However, when it acts 
on the helicities of the Dirac spinors as in the treatment 
by FKS, it is not the correct projector to ensure the spin- 
triplet states; the corresponding error modifies the Fermi 
motion effects of order (compare Eqs. (|41|l . 1421) with 
Eqs. (113-1201)). 

Now, we recognize that our vertex \2<j\i is different from 
that of RRML as well as that of FKS in two points: the 
correct spin structure with the projector TZ of Eq. 
and the off-shellness without the condition H39|l . Both 
of these points modify ImAl'^' of Eq. (|31|) at the order 
f ^ and lead to the additional Fermi motion corrections 
compared with those of RRML or FKS. 

There is, still, another "new" Fermi motion correction 
from a"^'^ of Eq. (HOJ): in almost all previous works in- 
cluding RRML and FKS, x appearing in 0-^^^ is replaced 
as 



(43) 



Because this x corresponds to the static limit of the on- 
energy-shell value of x (|18|) , the use of x with the explicit 
fc^- and a-dependence generates the O(f^) modification 
of lmM'-<\ 

Before ending this subsection, we mention the con- 
nection of our projector H22I) with the operator used in 
Ref. [231 to ensure the S'-wave cc states. When we im- 
pose the on-energy-shell condition 139() by hand, Eq. 126() 
reduces to 



2Af2^a(l-a) 



vxiq-kh^ei^^*S''''ux'{k), (44) 



with^^'' = g''" -2k''k'' /[mciMv+2mc)], usingthe Dirac 
equation for the spinors. This operator S'^'^ was used in 
Ref. 113. This result shows that Sf"" leads to the 0{v^) 
error, and it is not a projection operator, i.e., S^'^S^p ^ 



D. Production cross section with Fermi motion 
effects 

We are now in a position to present our final formula 
for the forward differential cross section of the diffractive 
heavy vector-meson production including all the Fermi 
motion effects discussed above, and also for the corre- 
sponding total production cross section. Substituting 
Eqs. (O, (1121, (EOl), I23-(|2n| into Eq. ^23, we get 
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eCcsQ 



Q — k 

:a(l - a)asiQlff)xG{x, Ql^^ 



k\ + mc{Mv + rric) 



a(l - a)M2 



as{Ql^)iG{x,Ql^) 



^ q;(1 - a) fc^ + Q 



(fci + Q 

Tic {{Mv + mjmc + a(l - a)M^} {k\ - Q^^ 



(45) 



2 { [a^ + (1 - af] Mv + m^} fclQ^ 



(46) 



r 



for the longitudinal and transverse polarizations, respec- 
tively. Here we have omitted the contributions that van- 
ish when integrating over the angle of fcj^ in Eq. 
The forward differential cross section is given as 



dcr(7(*)7V -> VN) 



dt 



1 



(C)|2 



(47) 



Following Refs. Q, we can calculate KeAi^''^ per- 

turbatively from lmM^''\ using ReA^'^'/ImA^'^' = p 
with 



^9 In {xG{x,Qls)) 



(48) 



The forward differential cross section H47|) is now given 
as 



da{-f^*^N VN) 



dt 



1 



^(ImTW 

C 



(C) 



(49) 

(The actual contribution from the real part is found to 
be less than 10% at HERA energies.) Combining the re- 
sult H49|) with the ^-dependence of the differential cross 
section observed in experiment, da /dt oc exp(i?v't), with 
a constant diffractive slope By 01 we can calculate the 
total cross section for the diffractive vector-meson pro- 
duction, 



cr(^(*)Ar ^ VN) 



1 dcr{j'^*^N 



VN) 



dt 



(50) 



To obtain ImTW^'^' of Eq. we use Eqs. (|45|l . I|46|) 

in Eq. Ij^ll) . In Sec. II C, we explained the Fermi motion 
effects by Taylor expanding Q^'^'>{a, k±) in the integrand. 
However, in the numerical calculation presented below, 
it is more convenient to evaluate directly the convolution 
integrals over a and k± without recourse to the Taylor 
expansion. The nonperturbative part of the light-cone 
WF, (f>*{a, k±), is constructed from a S-wave solution of 
the Schrodinger equation with a realistic potential be- 
tween c and c. We use the Cornell potential consisted of 
"Coulomb plus linear potential", V{r) = —k/t + r/a'^, 
with the parameters k and a chosen to reproduce the 
masses and the decay properties of the charmonia |l2j |. 
(The use of the other QCD-inspired potentials modifies 



the production cross section by at most 10%.) We follow 
the procedure of FKS in order to derive (l)*{a,k±) 
with the light-cone variables a, k± from the solution 
of the conventional "equal-time" Schrodinger equation, 
(f'NRilkl), where k denotes the (three-dimensional) vari- 
able for the space components of the momentum k: we 
employ the kinematical identification of the Sudakov 
variable a, which denotes the fraction of the "plus" com- 
ponent of the meson's momentum carried by the c-quark, 
with 



(51) 



This relation, together with the conservation of the over- 
all normalization of the WF, J dad'^k±\(j){a,k±)\^ /2 — 
J d'^k\(t)Nji{\k\)\'^ , allows us to express the light-cone WF 
in terms of the nonrelativistic WF as 





(52) 



The obtained WF (j>{a, k±) is in fact peaked at a = 1/2 
and k± = 0. When we take the static limit 0(a, fej^) 
const X(5(a- 1/2)5(2) and impose the condition H39|l , 
it is easy to see that the differential cross section (|49|) 
coincides with the old result derived by Ryskin 1], 



da 



48a, 



emQoff 



[asiQls)^Gix,Qls)y 



Ml 
(53) 

up to the small contribution due to p. Here, T^e stands 
for the decay width of the vector meson into an e^e~ pair. 
In the last parenthesis, the first and the second terms cor- 
respond to the production with transversely and longitu- 
dinally polarized photons, ctt and a^, respectively. As 
we will demonstrate in the next section, the dominant 
effects due to the Fermi motion corrections with the full 
WF H52|l modify these two contributions by the different 
factors that could in principle depend on as well as 
W. 
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On the other hand, if we formally suppose S> rr?^ in 
Eq. and impose the condition the result repro- 
duces the LLA formula of Brodsky tt al. |2| for the hard 

diffractive production. In this case, Q — ■ml+a{l—a)Q'^ 
becomes much larger than k]^ in Eqs. 1)451) . H46I) . except 
the end-point region of the a-integral, so that the a- 
dependence of the corresponding integrand plays a dom- 
inant role in the cross section compared with the k±- 
dependence; this situation corresponds to the light-cone 
dominance and the fej^-dependence is only responsible for 
the small higher-twist (~ corrections. However, in 

the region < to^ discussed below, the fei^-dependence 
is as important as the a-dependence in Eqs. H45|l . (|46|l . 
In particular, for the photoproduction — 0, the trans- 
verse quark motion represented by the fcj^-dependence 
could play a more dominant role in the Fermi motion 
effects. To show the quantitative role of the transverse 
quark motion, we will present the results of the cross sec- 
tions in the "static limit for the transverse motion" with 
the replacement in Eq. H31() . 

0(a,fei)^^(a),5(2)(fei), (54) 

where ip{a) = J (Pk±(j){a,k±), and compare with the 
results using the full WF . We will demonstrate that 
the significant effect from the transverse quark motion is 
a novel feature in the production of the heavy mesons. 



III. THE J/i/) AND V' PRODUCTION CROSS 
SECTION AND COMPARISON WITH DATA 

In this section, we present the numerical results, using 
our formulae (gHl) and ^SU^ with Eqs. (EU, (EHl), (gHIl, and 
(|52|1 . We show the cross section for the diffractive photo- 
and electroproductions of the heavy vector-mesons V = 
J/i^, t/j', and discuss the roles of the Fermi motion effects 
in detail. We use rric = 1.5 GeV for the charm quark 
mass, and Mj/^ = 3.10 GeV and M^, = 3.69 GeV for 
the mass of the vector mesons. For the slope parameter 
Bv of Eq. ISUJ), we adopt Bj/^ = 4.44 GeV^^ for J/^P 
and B^i = 4.31 GeV~^ for ip', that are extracted from 
the experiments l^lgHllf. 

For the gluon distribution function G{x, Q^ff) of 
Eqs. 601), we employ Gliick-Reya-Vogt (GRV) 

parametrization [l9j . Besides uncertainties in the knowl- 
edge of G{x,Q'^g), we cannot determine the scale Qlg, 
which enters into G{x, Q^g) as well as the strong cou- 
pling constant as((5^g), unambiguously within the accu- 
racy of the LLA. We also note that, within the LLA, the 
use of the NLO fits to G(i, Qlg) is not fully consistent. 
In fact, in the previous works, the various modifications 
of Ql fj were employed to try to sirnulate the important 
effects beyond the LLA HH 113, ill. In the present 
work, we do not pursue such "schemes" to fix the scale 
Qlff- In addition to the results with the "standard" scale 
= (Q2^A/^)/4 in Eqs. (gSJl and (gSJ), we will present 



the results with the replacement 

Qls ^ 2Ql^, (55) 

and thereby we demonstrate how much our results would 
be modified by the corresponding ambiguity. We note 
that actually the replacement (|55() gives very similar ef- 
fects to those of the "Q^ rescaling" employed by FKS, as 
one can read off from Fig. 9 in the first paper of Ref. 3]. 

A. Fermi motion effects and the gl/ut ratio 

As pointed out in Sec. II, the new Fermi motion effects 
in the LLA, compared with the previous works, originate 
from three points: 

(i) The spin WF for the ^5*1 cc state with the projector 

(ii) The off-shell cc ^ V vertex without using the con- 
dition 1(23). 

(iii) The modification of the gluon's longitudinal mo- 
mentum fraction probed by the process, as Eq. H18|l . 
due to the coupling with the cc pair in internal rel- 
ative motion. 

In order to understand the roles of these new effects, it 
is useful to compare our vertex involving (i) and (ii) with 
those of RRML and FKS: our vertex, Eqs. has 
quite different a and k± behavior from the RRML vertex 
(ESJ, lEH) as well as the FKS vertex glj, lO, for each 
(A, A') and for both longitudinally and transversely po- 
larized vector-mesons. Namely, (i) combined with (ii) is 
expected to produce strongly helicity-dependent Fermi 
motion effects. We further note that (ii) itself gives 
the helicity-dependent effects: by imposing the condition 
(|39|l on our vertex H26|) . we get Eq. 1441) . Apparently, 
this result corresponds to our vertex (|27|I - H3UI) with the 
term fc^ eliminated by using Eq. (|39|l . and this affects 

^-AA '-'^ ^'i- only. On the other hand, (iii) is not ex- 
pected to produce strong helicity-dependence: the effect 
of (iii) acts similarly on both longitudinal and transverse 
polarizations through the gluon distribution xG{x, Q'^g) 
in Eqs. (gHl) and 

To investigate the helicity-dependence of our new 
Fermi motion effects in detail, we compute the J /tp pro- 
duction cross section with longitudinally and transversely 
polarized photons (see Eqs. (j^ ). 

aL{l'^*^N ^ J/i) N) 
crj,(7(*)iV ^ J/lj} N) 




FIG. 3: The cross section for the diffractive J/tj} 
production with longitudinally polarized photon, as a 
function of at W = 97 GeV. The solid line shows 
our full result 15611 . including the new Fermi motion 
effects due to (i)-(iii). The dotted line is the result ex- 
cluding (iii), and the dashed line is the one excluding 
(ii) and (iii). The two-dot-dashed line shows the result 
without (i)-(iii), using the RRML vertex H35|l or the 
FKS vertex lUTl . 



FIG. 4: The cross section for the diffractive J/ip pro- 
duction with transversely polarized photon, as a func- 
tion of at W = 97 GeV. The solid line shows our 
full result 1571 , including the new Fermi motion effects 
due to (i)-(iii). The dotted line shows the result ex- 
cluding (iii), as well as excluding (ii) and (iii). The 
dot-dashed and the two-dot-dashed lines show the re- 
sult without (i)-(iii), using the RRML vertex H86^ and 
the FKS vertex 1421 . respectively. 



full result with (i)-(iii) _ ' 

without (iii) _ y ' 

. without (ii) and (iii) _ y ' 




(t [GeV^] 



FIG. 5: The ratio of the cross sections for the diffractive 
J/tp production with longitudinally and transversely polar- 
ized photons, as a function of at = 97 GeV. Each curve 
is obtained by taking the ratio of the corresponding results 
of Figs. El said |1] and the lines have the same meaning as in 
those figures. 



with Eqs. (gHl, and 



In Fig. 121 we show aL{'y*N — > J/i/; N) as a function 
of with the 7*-A^ center-of-mass energy fixed to be 
W = 97 GeV. The full result ^ is shown by the solid 
curve, which includes all the above-mentioned effects (i), 

(ii) , and (iii) . The dotted curve shows the result without 

(iii) , i.e., the result of Eq. (|5H|l using Eq. with the 
replacement ()4l-i(l . The dashed curve shows the result 
without (ii) and (iii). This is obtained from the result 
of the dotted curve by replacing the cc ^ V vertex by 
Eq. (|^ . By the two-dot-dashed curve, we also show the 
result without (i)-(iii): this is obtained from the result of 
the dashed curve by replacing the corresponding cc ^ V 
vertex (I44II by the RRML vertex (|35l) . which coincides 
with the FKS vertex for the longitudinal polarization (see 
Eq. l|lT|l ). The dotted curve is largely enhanced due to 
the effects of (i) and (ii), in comparison with the two- 
dot-dashed curve. Comparing the dotted and the solid 
curves, we see that (iii) gives the suppression. We note 
that the enhancement due to (ii) is almost canceled by 
(iii) in this channel (compare the solid and the dashed 
curves). Still, our full result shown by the solid curve is 
larger than the two-dot-dashed curve, demonstrating an 
important role of the Fermi motion effect by (i). 

The results for aT{j*N ^ J/tp N) are shown similarly 
in Fig. 12 as a function of Q"^ at W = 97 GeV. The full 
result H57|l involving (i)-(iii) is shown by the solid curve. 
The dotted curve shows the result with the replace- 
ment (|^ . corresponding to the case without (iii). We 
note that the dotted curve also corresponds to the case 
without (ii) and (iii), because our vertex f^ . (p?n|l coin- 
cide with the vertex I44|) for the transverse polarizations. 




FIG. 6: The total cross section for the electroproduc- 
tion, 'y* N ^ J/tp N, as a, function of at W = 90 
GeV, in comparison with the HERA data The 
sohd and the dotted hnes show our full result 1501 
with and without the replacement 11551 . respectively. 
The dot-dashed and the two-dot-dashed lines show the 
results without the effects due to (i)-(iii), using the 
RRML vertex ^ and the FKS vertex (gTJ, J^ . 
respectively. The dashed curve is obtained from the 
dotted curve with the replacement 11541 . 



FIG. 7: The total cross section for the photoproduc- 
tion, jN — > J/ip N, as a function of W . The lines 
have the same meaning as in Fig. Experimental 
data pmnts from the H10, ZEUS^^, E401i2§|, and 
E516|27ll experiments. The HERA data were measured 
in two leptonic decay channels, J/ip — > fi^ , and 
J/ip e^e~ . 



By the dot-dashed and the two-dot-dashed curves, we 
show the results without (i)-(iii): the dot-dashed (two- 
dot-dashed) curve is obtained from the result of the dot- 
ted curve by replacing the corresponding cc V vertex 
by the RRML vertex ^ (FKS vertex lO)- Compar- 
ing the sohd curve with the dotted curve, we see that (iii) 
gives the suppression, similarly to the case of (Ti of Fig. 13 
From the comparison of the dotted curve with the dot- 
dashed curve, we may say that the effect due to (i) gives 
rise to the suppression, contrary to the case of ctl; recall 
that the RRML vertex does not take into account the 
projector TZ for the 5- wave cc state at all, i.e., TZ ^ \. 
We note that the two-dot-dashed curve using the FKS 
vertex is much more suppressed than the dotted curve. 
This indicates that the FKS projector H4U|) . which has 
an oversimplified structure omitting the A' = — A compo- 
nent of the vertex as in Eq. H42|l. overestimates the effect 
of (i) due to the spin WF. 

To summarize the helicity-dependence of the Fermi 
motion effects due to (i)-(iii): the effect due to (i) en- 
hances (Tl, while it suppresses ct; the effect due to (ii) 
enhances ctl, and it does not change ctt; the effect due to 
(iii) suppresses both ctl and or. We also note that the 
"oversimplified" FKS projector l|40|l significantly overes- 
timates the suppression of ctt due to (i) . 

Those strong helicity-dependence from (i), (ii) gives 
rise to notable behavior of the ratio (7l(7*A^ 
J/^ N)/aT{l*N J/i! N), which is shown in Fig. 
Here each curve is obtained by taking the ratio of the 
results expressed by the corresponding curves in Figs. |2| 



and01 Note that, if we plotted the dashed curve in Fig.^ 
(the dot-dashed curve in Fig. |2J), it would completely 
coincide with the dotted (two-dot-dashed) curve. The 
dashed curve is enhanced in comparison with the dot- 
dashed curve due to the "combined" effect of (i) on the 
ratio (iL/aT, and it is further raised to the dotted curve 
due to the effect of (ii). Comparing the dotted and the 
solid curve, we see that (iii) gives the enhancement, but it 
is not significant; the effects of (iii) affect ctl and ctt sim- 
ilarly, so that they largely cancel for glIot- Note the 
quite strong enhancement of the two-dot-dashed curve 
compared with the dot-dashed curve, because of the sig- 
nificant overestimate of the suppression of gt due to (i) . 

Another interesting aspect of Figs. lSISl is the character- 
istic behavior as a function of . Rather different shape 
between ct^ and ut of Figs. 13 and 0] corresponds to the 
behavior of the ratio (JlIut in Fig. [3 which increases al- 
most linearly with increasing . As is well-known, this 
Q^-dependence is due to the behavior of the light-cone 
WF for the longitudinal photon, Eq. (|12() . whose first 
term is proportional to Q (see also Eq. (|45|) '). Also, if we 
compare the explicit a- or fe_L-dependent terms between 
Eqs. H45() and (|46() . those terms reflecting (i), (ii) depend 

2 

differently on Q . As we have shown above, however, 
the Fermi motion effects, produced via the convolution 
of Eq. (PJ) using Eqs. do not significantly 

modify the shape of ctl, ctt, and ul/ut- The Fermi mo- 
tion effects mostly give overall enhancement or suppres- 
sion, by the different factors for gl and gt respectively. 
In the next subsection, the ratio gl/gt, as well as the 
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corresponding total cross section cr = ctl + o^t, will be 
compared with the recent HERA data. 



B. The J/ip photoproduction and electroproduction 
cross section 

Now, we present our numerical results for the J/^p 
photo- and electroproduction cross section (t{j'-*'>N 
J/ip N) (see Eq. (|5U|) '). including the new Fermi motion 
effects due to (i)-(iii), and make a comparison with the 
available data. 

In Figs. El and we show a{-f*N J/^ N) as a 
function of ^ith = 90 GeV fixed, and uinN 
J I N) for the photoproduction (Q^ = 0) as a function of 
W, respectively. In both figures, the dotted curve shows 
our result (EOJ with Eqs. (EU, (03, llHI), and (|S2Il- The 
solid curve is obtained by applying the replacement H55|) 
to the dotted curve. The dot-dashed and two-dot-dashed 
curves show the results without the effects due to (i)-(iii), 
using the RRML vertex (jSHj), JSEI) and the FKS vertex 
()41|l . (|42|l for the effective cc ^ V vertex, respectively. 
Note that the dotted, dot-dashed, and two-dot-dashed 
curves correspond to the cases of the solid, dot-dashed, 
and two-dot-dashed curves in Figs. 0113 respectively. We 
also show, by the dashed curve, the result obtained from 
the dotted curve with the replacement H54|) . in order to 
demonstrate the role of the Fermi motion effects in the 
transverse direction. 

For < 20 GeV'^, ctt is dominant in comparison 
with aL (see Figs. \'6l6\i . Therefore, the behavior of the 
curves of Fig. El in this region is very similar to that of 
Fig. ^ The behavior due to <tl shows up for > 20 
GeV^, raising the tail. In particular, such characteristic 
behavior is clearly seen for the solid, dotted, and dot- 
dashed curves in Fig. EJ which are in agreement with 
the data. We note that the corresponding three curves 
reproduce the behavior of the data also in Fig. [7\ 

One observes that the curves in Figs. El El have some 
common features: 

(a) For most region of and W shown in those fig- 
ures, our full result (dotted curve) shows slight (10 ~ 
20%) suppression compared to the dot-dashed curve, i.e., 
the case without (i)-(iii). This reflects the correspond- 
ing suppression in Fig. 0] because cr « ctt for < 20 
GeV^ in Fig. El and ct = or for Fig.Q We note that the 
dot-dashed curve of Fig. [7| coincides with the estimate of 
the Fermi motion effects for the J/tp photoproduction by 
RRML 4], up to the different form of the nonperturba- 
tive part of the hght-cone WF, (/)*(«, fc^). 

(b) The two-dot-dashed curve is suppressed very 
strongly in comparison with the other curves. This again 
reflects the corresponding strong suppression in Fig. 01 
due to the "overestimate" for the effect of (i) by using the 
FKS projector (gnj. We note that the two-dot-dashed 
curve in Figs. El U\ coincides with the result calculated 
by FKS up to some elaborate corrections considered 



in their paper, such as the "Q^ rescaling", the running 
quark mass, etc. 

(c) The replacement (|55f) pushes up the dotted line to 
the solid line, by a factor of 1.1 ~ 1.5 for Fig. El and 
1.2 ~ 1.8 for Fig. El One may say that the ambiguity for 
the scale Q^g would not lead to so significant modifica- 
tion of the results for W < 100 GeV. For higher W, the 
corresponding enhancement of the results becomes pro- 
nounced, because the scale-dependence of the gluon dis- 
tribution G{x, Qlg) becomes stronger for smaller x (see 
Eq. (|18|l ). We note that, if one applies the replacement 

to the dot-dashed or the two-dot-dashed curves, 
those curves are enhanced by almost the same factor as 
the dotted curve. Thus, the two-dot-dashed curve could 
be consistent with the data, too. This might suggest that 
the NLO perturbative corrections could give the effects 
of the same order as the Fermi motion corrections to the 
cross section for the charmonium production. 

(d) Comparing the dotted and the dashed curves, we 
recognize that the total contribution of the transverse 
Fermi motion gives very strong suppression by a factor 
~ 1/4. This clearly shows the significant role of the trans- 
verse motion of the quarks for the charmonium produc- 
tion (see the discussion above Eq. ifS^ I. 

Fig. El shows the J/^ electroproduction cross section, 
(t(7*7V ^ J/V' N), as a function of W for = 3.5^ iq.i^ 
and 33.6 GeV^. The solid and the dotted lines have the 
same meaning as in Fig. (3 The W dependence of the 
results is similar to the photoproduction case in Fig. El 
and both solid and dotted curves reproduce the data. 
Although we have not shown explicitly, we have also cal- 
culated the results corresponding to the dot-dashed as 
well as the two-dot-dashed curve of Fig. El and find that 
those results are also consistent with the data in Fig. El 

Fig. El compares the results of (Tl/ct of Fig. El with 
the recent HERA data for = 97 GeV |1 1|. The 
solid, dash-dotted, and two-dot-dashed curves are iden- 
tical with the corresponding curves in Fig. El respec- 
tively. The dashed curve is obtained by applying the 
replacement (|54f) to the solid curve, while the dotted 
line shows the result corresponding to Eq. H53|) . i.e., 
fL/o'T= Q^/Mj^^. Thus, the total contribution due to 
the Fermi motion is significant, and the quantitative role 
of the transverse quark motion is similar to that of the 
longitudinal one. 

The results with the replacement (|55|) are not shown 
in Fig. El We actually find that, for the ratio (Jl/(Jt, 
the effects due to the replacement (|55|) almost cancel be- 
tween the numerator and the denominator. This suggests 
that the ratio (Tl/ctt would allow a "clean" theoretical 
prediction insensitive to the ambiguity for the scale Q^g- 
One might further expect that (Tl/'^t could be insen- 
sitive even to other ambiguities or corrections like the 
uncertainties in the i-dependence of G{x, Q^g) and the 
NLO perturbative corrections to the diffractive ampli- 
tude. Moreover, as demonstrated in Sec. HI A, the strong 
helicity-dependence of the new Fermi motion effects due 




W [GeV] Q' [GeV^] 



FIG. 8: The total cross section for the electropro- 
duction, y'N J/ij) N, as a function of W for 
= 3.5,10.1, and 33.6 GeV^, in comparison with 
the recent HERA data 0,0. The sohd and the dotted 
hnes have the same meaning as in Fig. |7| The open 
squares show the ZEUS data at = 13 GeV^, which 
is somewhat larger than the indicated value of . 



FIG. 9: The ratio (Ti^/gt for the diffractive J/ip pro- 
duction in comparison with the HERA data 5, 6] . The 
solid, dash-dotted, and two-dot-dashed curves are iden- 
tical with the corresponding curves in Fig. |K| respec- 
tively. The dashed curve is obtained by applying the 
replacement 1541 to the solid curve, while the dotted 
line shows the result corresponding to Eq. (1531 . 



to (i), (ii) modifies (JlIut strikingly, so that we expect 
that the ratio cr^ /ctt could be the most suitable discrim- 
inator among the predictions of the Fermi motion effects 
by various models. Unfortunately, however, the available 
data in Fig. [U] are not enough for this purpose. 

C. The ?/;' photoproduction and electroproduction 
cross section 

It is straightforward to extend our calculation in 
Sec. Ill B to the i/^' production, by substituting the 
solution for the first radial excitation of the relevant 
Schrodinger equation into (/>jv_R of Eq. H52|) . and by the 
trivial replacements, Mj/^ M^i, Bj/^ — > -B^', in 
the relevant formulae. We discuss our numerical re- 
sults for the V'' photo- and electroproducion cross sec- 
tion a [^i*)N ip'N), including the new Fermi motion 
effects due to (i)-(iii), and also make a comparison with 
the available data on the ratio of ij}' to J /"ip cross section. 

In Fig. ^1 we show the ratio a{'^*N 
'tp'N)/a{j*N — + J/ip N) for the electroproduction as a 
function of at W = 90 GeV. The lines have the same 
meaning as in Fig. |SJ the solid curve shows the ratio of 
a{'~f*N ^'N), which is calculated similarly to the solid 
curve in Fig. El to a{'~f*N J/ip N), which is given by 
the solid curve of Fig. El the dotted curve shows the ratio 
of a{j*N ip'N), which is calculated similarly to the 
dotted curve in Fig. O to a{j*N — > J/ip TV), which is 
given by the dotted curve of Fig. and so on. 

We see that the solid curve almost coincides with the 
dotted curve, because, similarly to the ratio gl/i^t of 
Fig. El the effects due to the replacement (|55|l almost can- 



cel between the numerator and the denominator. Note 
that, for the dotted curve of the present case, we use 
Ql^ = (g2 + Mj/^)/4 for a,{Ql^)G{iMl^) of the de- 
nominator, and Q^^j = (Q^ -|- M|,)/4 for that of the nu- 
merator (see Eqs. (|45ll . (|46|l ^: actually, this difference of 
Qcff between J/ -0 and gives the small difference of the 
solid and the dotted curves for the low region. Ap- 
parently, such small difference is irrelevant, and for high 
Q^, i.e., ^ Mj^^,M^,, the difference between those 
two curves disappears. However, the corresponding dif- 
ference of Q^ff could produce relevant behavior for other 
observables (see the discussion about Fig. Illlbelowl. 

The comparison of the dashed line with the dotted 
line suggests that the Fermi motion effects are more pro- 
nounced for the ijj' production than for the J/ip one. In 
fact, the typical velocity v of Eq. H33|) is larger for the 
"25-state" ^' than for the "l^-state" J/i). We also rec- 
ognize an interesting point: only the dashed line has the 
"flat" behavior, while all the other lines have qualita- 
tively a similar shape with each other, decreasing steeply 
with decreasing Q^. It is easy to see that this latter be- 
havior is characteristic of the Fermi motion effects for 
the case of the radial excitation, reflecting that the cor- 
responding meson WF has a "node" : for <T{'y*N — > iIj'N), 
we calculate the convolution H31() with the nonperturba- 
tive WF (/)*(a,fc^) for ^' , which is given as Eq. in 
terms of the nonrelativistic, 25-state WF (pNR- Because 
(f>NR. has the node at r ^ l/(mcw) in the coordinate space, 
the convolution 1)311) suffers from the strong cancellation. 
We note that the behavior of Q^''\a,k±) as a function 
of a,k± is controlled by the scale (see Eqs. (j^ . 
(|46ll ): therefore, the cancellation can be avoided when 
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FIG. 10: Same as Fig. |S|but for the ratio of the total 
cross section for the electroproduction, 7* A/' — > V'-^i 
to that for 7*iV ^ J/ip N. 



FIG. 11: Same as Fig. |7|but for the ratio of the total 
cross section for the photoproduction, 7A'^ — > ip' N , to 
that for 7A'^ ^ J/ip N. Ex peri mental data points from 
earlier fixed target |33 IH 113 and the HI US 

ITU experiments. 



3> m^, i.e., Q^^") is "localized" sufficiently in the co- 
ordinate space representation. This effect leads to the 
steep Q^-dependence for a{'y*N tp'N), which is ab- 
sent for a{j*N "> J/ip N), so that one obtains the be- 
havior of the curves in Fig. ^1 The reason why the 
dashed curve has the flat behavior is that we use the 
replacement (|54|) for the ip' WF in the numerator, as 
well as for the J/4> WF in the denominator, to calculate 
cr(7*7V i}'N)/a{-/*N -> J/^j N) in this case. 

Our full result (dotted or solid line), including the new 
Fermi motion effects due to (i)-(iii), is suppressed by a 
factor of about 2/3 than the dot-dashed line, which cor- 
responds to the case without (i)-(iii) and uses the RRML 
vertex (|35|l . On the other hand, the two-dot-dashed 

line, which corresponds to another case without (i)-(iii), 
shows much stronger suppression. This is due to the com- 
bined effect of the "overestimate" of the effect of (i) by 
using the FKS projector (|40|l . and of the fact that the 
Fermi motion effects are more pronounced for the radial 
excitation ip' than for J/^. 

In Fig.ini we show the ratio ai-jN i)'N)/a{^N -> 
J/ip N) for the photoproduction (Q^ = 0) as a function 
of W. The lines have the same meaning as in Fig.^] For 
almost the whole region of W, we see the behavior simi- 
lar to that observed for ^ in Fig.^l In particular, 
the dashed curve with the replacement (|54|l largely over- 



estimates the data by a factor of 2 '--^ 3, in contrast with 
the other curves. This suggests that the strong trans- 
verse Fermi motion effects, reflecting the radial shape of 
the 25'-state WF, is essential to explain the tendency of 
the data. On the other hand, the two-dot-dashed curve, 
using the FKS vertex H41|) , (|42|l , underestimates the data 
by a factor of 1/3 ~ 1/2. We will add some comments 
on this point in the final part of this subsection. 

Our full result, as well as all the other curves in Fig. II II 
shows a slight increase with increasing W, above ~ 10 
GeV, and this behavior seems to match with the ten- 
dency of the recent HI data the present cal- 
culation, this weak W dependence originates from the 
difference of Q^^j = {Q^ + My) /A in the gluon distribu- 
tion G(i,(3pff), corresponding to F = -0' and J/i/': the 
larger Q^jj gives the faster increase of G{x, Q^g) with de- 
creasing X (see Eq. ()18(l 'l. Actually, this fact motivates 
our choice as Q^g = {Q^ + My) /A following Ref. al- 
though, within the LLA, other choice is possible. For 
example, in Ref. Q^g — {Q^ + 4m^)/4 was used. 
We have checked that this choice makes the W depen- 
dence of the ratio a{jN ^p'N)/a{jN J/i/' N) al- 
most constant. Another choice used in the literature is 
Qlfi = a(l — a)Q^ + ml + k^j_ 0|. In this case, we find 
that the ratio aijN ip' N)/(j{'yN J/ip N) shows a 
weak rise similarly to that in Fig. ^2 due to the convo- 
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FIG. 12: Same as Fig. |H|but for the ratio of the total FIG. 13: Same as Fig. |n|but for the diffractive xp' pro- 

cross section for the electroproduction, 7* A/' — > ^p'N, duction. 
to that for 7*iV ^ J/ip N. 



lution with the ip' or J/ip WF 4>*{a, k±) with respect to 

A similar rise as a function of W is observed in Fig. 1121 
which shows the ratio a{j*N tp'N)/a{-f*N J/ip N) 
for the electroproduction for = 3.5,10.1, and 33.6 
GeV^. The solid and the dotted lines have the same 
meaning as those in Fig. ITTI 

Now, we go over to the helicity dependence in the 
diffractive production. Fig. El shows the ratio 
(Tl{i*N ip'N)/aT{l*N ^'N) for the electropro- 
duction of %p' as a function of Q'^ dX W = 97 GeV. The 
lines have the same meaning as in Fig. |51 for the J /tp 
production. We see that the strong helicity-dependence 
observed in Fig. |51 is even more pronounced for the -0' 
production. We expect that such behavior is insensitive 
to elaborate corrections, and ctl/ct for the ip' produc- 
tion allows a clean prediction, although we do not have 
the experimental data to be compared. 

We also show, in Fig. the W dependence of the 
ratio (Jl{i*N ^/j'N)/aT{j*N ip'N) by the dotted 
lines, and of ctl(7*7V -> J/i) N)/(Jt{i*N Jl^ N) by 
the solid lines, for = 10,20, and 90 GeV^. All lines 
are calculated using our full results (|^ . (|57|l with the 
replacement . We see that those results are almost in- 
dependent of W for the energy region available at HERA. 
This implies that the relevant W dependence through x 
of Eqs. (|45|l . H4()|) cancels with each other between the 
denominator and the numerator of (Tl/ctt, even though 



we take into account the a and k± dependence of x as 
in Eq. H18|) . When we make the replacement H43|) . the 
approximately constant behavior of the /ct ratio as a 
function of W would be an immediate consequence of the 
cancellation of the gluon distribution G{x, Q^ff ) between 
the denominator and the numerator, as is also seen from 
the corresponding result of Ref. [s^ ■ 

Finally, we mention some recent works which also dis- 
cussed the ratio of ip' to J/tp cross sections. In Figs. 1101 
111! we have shown the strong suppression of the two- 
dot-dashed curve compared with the data. Such strong 
suppression with the use of the FKS vertex H41|l . 142|) 
was first claimed for the photoproductions in Ref. (2Sj| . 
In the present work, we have clarified that the main rea- 
son of this too strong suppression is the "overestimate" 
of the Fermi motion effect of (i) by using the "oversimpli- 
fied" FKS projector H40|) . So far, there exist some works 
[2^ 01 , which suggested other mechanism to avoid the 
corresponding strong suppression. Ref. considered 
the role of the "new" spin structure for the cc state 
originating from the IVIelosh spin rotation, which is to 
relate the quantities in the constituent quark model to 
those in the infinite momentum frame. Ref. dis- 
cussed the role of the higher order terms in the power 
series of • b in the integrand of Eq. 1)16(1 . which corre- 
sponds to going beyond the LLA, i.e., the color-dipole 
picture. The detailed comparison of our results with 
those of Refs. [23ll23| is beyond the scope of this work. 
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FIG. 14: The ratio of the cross sections for the diffractive 
vector-meson production with longitudinally and transversely 
polarized photons, as a function of W at — 10, 20, and 90 
GeV^. The solid lines show the J/ip production, while the 
dotted lines show the tp' production. All lines are calculated 
using Eqs. 15611 . 1571 with the replacement 1551 . 



IV. SUMMARY AND DISCUSSION 

In the present paper, we have reexamined the Fermi 
motion corrections to the diffractive photo- and electro- 
productions of the heavy vector-mesons, J/ip and ip', in 
the LLA of pQCD. We have taken into account all the 
Fermi motion corrections arising from the relative mo- 
tion of quarks inside the charmonium, which is treated 
as a nonrelativistic bound state of c and c. The key in- 
gredients in our approach are the projector TZ to ensure 
the spin structure for the cc bound state, the off- 
shellness in the hadronization (cc V) vertex, and the 
modification of the Bjorken-x, i.e., the gluon's longitu- 
dinal momentum fraction probed by the process, due to 
the coupling with the cc pair in internal relative motion. 
Although these three contributions were not considered 
properly in the previous works, we have demonstrated 
that all these contributions should be included in the 
LLA diffractive amplitude in QCD, and that all of them 
produce the new Fermi motion effects of order with v 
being the heavy-quark velocity inside the vector meson. 

We have presented a quantitative estimate of the Fermi 
motion effects on the diffractive production cross sec- 
tions, including our three new contributions. In the QCD 
factorization formula for the diffractive charmonium pro- 
duction, the nonperturbative, internal motion of quarks 
is represented by the light-cone WFs for the charmonium, 
which we have constructed from the corresponding S- 
wave solution of the Schrodinger equation with a realis- 
tic potential between c and c. We find that the net effect 



from the Fermi motion gives moderate suppression of the 
total cross section for the J/ ip photo- and electroproduc- 
tions. The corresponding suppression is similar to that 
obtained by RRML p| for the J/^p photoproduction, but 
it is much weaker than that observed by FKS 3] for the 
photoproduction as well as the electroproduction. We 
have also clarified the strong helicity-dependence of our 
new Fermi motion effects, and predicted novel behavior 
in the helicity-dependent cross sections: the slope of the 
longitudinal to transverse production ratio (Tl/<tt as a 
function of distinguishes clearly our result from the 
results corresponding to the previous works. Because the 
'^l/o't ratio is expected to be insensitive to elaborate cor- 
rections or theoretical uncertainties, we propose it as a 
clean discriminator among the various predictions of the 
Fermi motion effects. For the detailed comparison of the 
(^l/ct ratio between theory and experiment, we should 
wait for the new precise data, especially for the high 
region. 

We have also observed the similar behavior for the ip' 
photo- and electroproductions. Our results for the ratio 
of to J/ip cross sections indicate that the suppression 
of the cross section becomes somewhat stronger for the 
ip' production than for the J/ip one, reflecting rapider 
motion of quarks inside the first radial excitation. We 
have demonstrated that the ip' to J/ tp ratio is also useful 
in discriminating the various theoretical calculations of 
the Fermi motion effects. Another interesting point is 
that the dependence of the ip' to J/ ip ratio is sensitive 
to the detailed shape of the WFs of the ip' as well as of 
the J/ip, while it is insensitive to the behavior of the 
gluon distribution. 

As a novel feature in the heavy meson production, we 
have emphasized that the Fermi motion in the transverse 
direction is as important as that in the longitudinal di- 
rection. The transverse quark motion plays a major role 
to produce characteristic effects due to the Fermi motion, 
especially the strongly helicity-dependent effects. 

The comparison of our results with the HERA data 
has been made in this paper. Taking into account the 
theoretical ambiguity associated with the LLA, we may 
conclude that our full results involving the new Fermi 
motion effects are in agreement with the recent data. 
However, it should be taken with reservation, because, in 
this work, we have not pursued some corrections, which 
could in principle produce the effects of the same order 
as the Fermi motion corrections. Those include genuine 
relativistic corrections due to the "small components" of 
the meson WF , and the next-to-leading perturbative 
corrections to the photon WF ^I'''' and the cc-gluon hard 
scattering amplitude yl'^'^f (see Eqs. (Q, Q). Leaving 
aside those unresolved problems for the present, the com- 
parison has shown that our results reproduce the overall 
behavior of the data as functions of as well as the 
energy W, over a range of observables. Therefore, what 
we can learn from this study is that the Fermi motion ef- 
fects provide natural mechanisms within the LLA, which 
give the characteristic behavior of the total cross sec- 
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tions as well as of the helicity-dependent cross sections 
observed in experiment. For the more detailed quanti- 
tative comparison with the data, one would eventually 
have to take into account further sophisticated correc- 
tions besides above-mentioned ones: those include, e.g., 
the contribution of higher Fock states such as \ccg) in 
the charmonia, the "skewedness" of the gluon distribu- 
tion 34], and the VF-dependence of the t-slope parameter 
By for the differential cross section da/dt p35l |. 
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